(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     11314,        355]
NotebookOptionsPosition[     10316,        315]
NotebookOutlinePosition[     10669,        331]
CellTagsIndexPosition[     10626,        328]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}], "\[IndentingNewLine]", 
 RowBox[{"a", "=", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"m1", " ", "v10"}], "+", 
        RowBox[{"m2", " ", "v20"}]}], "\[Equal]", 
       RowBox[{
        RowBox[{"m1", " ", "v1"}], "+", 
        RowBox[{"m2", " ", "v2"}]}]}], ",", 
      RowBox[{
       RowBox[{
        RowBox[{"m1", " ", 
         RowBox[{"v10", "^", "2"}]}], "+", 
        RowBox[{"m2", " ", 
         RowBox[{"v20", "^", "2"}]}]}], "\[Equal]", 
       RowBox[{
        RowBox[{"m1", " ", 
         RowBox[{"v1", "^", "2"}]}], "+", 
        RowBox[{"m2", " ", 
         RowBox[{"v2", "^", "2"}]}]}]}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"v1", ",", "v2"}], "}"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.9490449446317883`*^9, 3.949044952949711*^9}, {
  3.9490449891714177`*^9, 3.9490450604540195`*^9}, {3.9490451335365725`*^9, 
  3.9490451532152586`*^9}, {3.9490471949684787`*^9, 3.949047195613185*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"v1", "\[Rule]", "v10"}], ",", 
     RowBox[{"v2", "\[Rule]", "v20"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"v1", "\[Rule]", 
      FractionBox[
       RowBox[{
        RowBox[{"m1", " ", "v10"}], "-", 
        RowBox[{"m2", " ", "v10"}], "+", 
        RowBox[{"2", " ", "m2", " ", "v20"}]}], 
       RowBox[{"m1", "+", "m2"}]]}], ",", 
     RowBox[{"v2", "\[Rule]", 
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "m1", " ", "v10"}], "-", 
        RowBox[{"m1", " ", "v20"}], "+", 
        RowBox[{"m2", " ", "v20"}]}], 
       RowBox[{"m1", "+", "m2"}]]}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.9490471969185257`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"v1", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"a", "[", 
     RowBox[{"[", "2", "]"}], "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}], "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"v2", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"a", "[", 
     RowBox[{"[", "2", "]"}], "]"}], "[", 
    RowBox[{"[", "2", "]"}], "]"}], "[", 
   RowBox[{"[", "2", "]"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.9490450788721924`*^9, 3.949045097342475*^9}, {
  3.949045212930853*^9, 3.9490452897163153`*^9}, {3.949045338359805*^9, 
  3.9490453387745275`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"m1", " ", "v10"}], "-", 
   RowBox[{"m2", " ", "v10"}], "+", 
   RowBox[{"2", " ", "m2", " ", "v20"}]}], 
  RowBox[{"m1", "+", "m2"}]]], "Output",
 CellChangeTimes->{
  3.949045343627777*^9, 3.9490455464442263`*^9, 3.949046888459261*^9, 
   3.9490470410815926`*^9, {3.9490471477579207`*^9, 3.9490471666176567`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"2", " ", "m1", " ", "v10"}], "-", 
   RowBox[{"m1", " ", "v20"}], "+", 
   RowBox[{"m2", " ", "v20"}]}], 
  RowBox[{"m1", "+", "m2"}]]], "Output",
 CellChangeTimes->{
  3.949045343627777*^9, 3.9490455464442263`*^9, 3.949046888459261*^9, 
   3.9490470410815926`*^9, {3.9490471477579207`*^9, 3.9490471666226597`*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"m2", "=", 
     RowBox[{"m", " ", "m1"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"v20", "=", 
     RowBox[{"n", " ", "v10"}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.9490452548531528`*^9, 3.9490452612931643`*^9}, {
   3.9490453640650873`*^9, 3.9490453862126675`*^9}, 3.9490455207206583`*^9, {
   3.949046803378358*^9, 3.9490468069537067`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"v1", "=", 
   RowBox[{"v1", "//", "Simplify"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"v2", "=", 
  RowBox[{"v2", "//", "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n", "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Plot", "[", 
  RowBox[{"v2", ",", 
   RowBox[{"{", 
    RowBox[{"m", ",", "0.1", ",", "2"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.9490453974199243`*^9, 3.949045401855633*^9}, {
  3.9490454983756237`*^9, 3.9490455084353695`*^9}, {3.9490468186162148`*^9, 
  3.9490468722926674`*^9}, {3.949047092982633*^9, 3.949047097471572*^9}, {
  3.949047156784661*^9, 3.9490471570715237`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"2", " ", "m1", " ", "v10"}], "-", 
   RowBox[{"m1", " ", "v20"}], "+", 
   RowBox[{"m2", " ", "v20"}]}], 
  RowBox[{"m1", "+", "m2"}]]], "Output",
 CellChangeTimes->{{3.949045509228464*^9, 3.9490455560104456`*^9}, {
   3.9490468277774696`*^9, 3.9490468737234955`*^9}, 3.94904704690952*^9, 
   3.9490470988909187`*^9, {3.9490471581833534`*^9, 3.9490471706380744`*^9}}],

Cell[BoxData[
 GraphicsBox[{},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{0.5, 0},
  Method->{},
  PlotRange->{{0.1, 2}, {0., 0.}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{{3.949045509228464*^9, 3.9490455560104456`*^9}, {
   3.9490468277774696`*^9, 3.9490468737234955`*^9}, 3.94904704690952*^9, 
   3.9490470988909187`*^9, {3.9490471581833534`*^9, 3.9490471706420755`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"DiscretePlot3D", "[", 
    RowBox[{
     FractionBox[
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "m"}], ")"}], " ", "n"}]}], ")"}], " ", 
       "v10"}], 
      RowBox[{"1", "+", "m"}]], ",", 
     RowBox[{"{", 
      RowBox[{"m", ",", "0", ",", "2"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "0", ",", "0"}], "}"}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"params", ",", 
     RowBox[{"-", "2"}], ",", "2"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.9490467029586897`*^9, 3.949046786166475*^9}},
 NumberMarks->False],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`params$$ = -2, Typeset`show$$ = True, 
    Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", 
    Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = 
    "\"\:65e0\:6807\:9898\"", Typeset`specs$$ = {{
      Hold[$CellContext`params$$], -2, 2}}, Typeset`size$$ = {
    360., {142., 146.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True, $CellContext`params$1012$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, "Variables" :> {$CellContext`params$$ = -2}, 
      "ControllerVariables" :> {
        Hold[$CellContext`params$$, $CellContext`params$1012$$, 0]}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> 
      DiscretePlot3D[(
         2 + (-1 + $CellContext`m) $CellContext`n) ($CellContext`v10/(
         1 + $CellContext`m)), {$CellContext`m, 0, 2}, {$CellContext`n, 0, 
         0}], "Specifications" :> {{$CellContext`params$$, -2, 2}}, 
      "Options" :> {}, "DefaultOptions" :> {}],
     ImageSizeCache->{411., {188., 193.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{
  3.9490455927089725`*^9, {3.9490467334428263`*^9, 3.9490467873568473`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["v2"], "Input",
 CellChangeTimes->{{3.9490470650971713`*^9, 3.9490470654993534`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"2", " ", "v10"}], 
  RowBox[{"1", "+", "m"}]]], "Output",
 CellChangeTimes->{3.9490470662570105`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["a"], "Input",
 CellChangeTimes->{3.949047177696017*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{"m1", " ", "v10"}], "-", 
        RowBox[{"m2", " ", "v10"}], "+", 
        RowBox[{"2", " ", "m2", " ", "v20"}]}], 
       RowBox[{"m1", "+", "m2"}]], "\[Rule]", "v10"}], ",", 
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "m1", " ", "v10"}], "-", 
        RowBox[{"m1", " ", "v20"}], "+", 
        RowBox[{"m2", " ", "v20"}]}], 
       RowBox[{"m1", "+", "m2"}]], "\[Rule]", "v20"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{"m1", " ", "v10"}], "-", 
        RowBox[{"m2", " ", "v10"}], "+", 
        RowBox[{"2", " ", "m2", " ", "v20"}]}], 
       RowBox[{"m1", "+", "m2"}]], "\[Rule]", 
      FractionBox[
       RowBox[{
        RowBox[{"m1", " ", "v10"}], "-", 
        RowBox[{"m2", " ", "v10"}], "+", 
        RowBox[{"2", " ", "m2", " ", "v20"}]}], 
       RowBox[{"m1", "+", "m2"}]]}], ",", 
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "m1", " ", "v10"}], "-", 
        RowBox[{"m1", " ", "v20"}], "+", 
        RowBox[{"m2", " ", "v20"}]}], 
       RowBox[{"m1", "+", "m2"}]], "\[Rule]", 
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "m1", " ", "v10"}], "-", 
        RowBox[{"m1", " ", "v20"}], "+", 
        RowBox[{"m2", " ", "v20"}]}], 
       RowBox[{"m1", "+", "m2"}]]}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.949047178882608*^9}]
}, Open  ]]
},
WindowSize->{707, 504},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (2013\:5e742\:670813\
\:65e5)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 1059, 29, 92, "Input"],
Cell[1641, 53, 749, 23, 77, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2427, 81, 597, 17, 52, "Input"],
Cell[3027, 100, 366, 9, 48, "Output"],
Cell[3396, 111, 366, 9, 48, "Output"]
}, Open  ]],
Cell[3777, 123, 449, 11, 72, "Input"],
Cell[CellGroupData[{
Cell[4251, 138, 668, 15, 92, "Input"],
Cell[4922, 155, 418, 9, 48, "Output"],
Cell[5343, 166, 525, 14, 241, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5905, 185, 729, 23, 75, "Input"],
Cell[6637, 210, 1737, 35, 398, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8411, 250, 98, 1, 31, "Input"],
Cell[8512, 253, 140, 4, 48, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8689, 262, 69, 1, 31, "Input"],
Cell[8761, 265, 1539, 47, 196, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
